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Collisionless, cold dark matter in the form of weakly- interacting massive particles (WIMPs) is well- 
motivated in particle physics, naturally yields the observed relic density, and successfully explains 
structure formation on large scales. On small scales, however, it predicts too much power, leading to 
cuspy halos, dense cores, and large numbers of subhalos, in apparent conflict with observations. We 
consider superWIMP dark matter, produced with large velocity in late decays at times 10 s s — 10 s s. 
As analyzed by Kaplinghat in a more general setting, we find that superWIMPs have sufficiently 
large free-streaming lengths and low phase space densities to help resolve small scale structure 
problems while preserving all of the above-mentioned WIMP virtues. 

PACS numbers: 95.35.+d, 95.30.Cq, 98.80.Cq, 12.60.-i 



The microscopic identity of dark matter (DM) is one 
of the major puzzles in basic science today. In the cur- 
rent standard cosmological picture, the Universe contains 
non-baryonic dark matter with abundance ^dm^ 2 = 
0.095 — 0.129 1], where Hdm is the energy density in units 
of the critical density, and h ~ 0.71 is the reduced Hubble 
parameter. This component is typically assumed to be 
cold, collisionless, and non-self-interacting dark matter, 
which we refer to as CDM throughout this work. CDM 
is remarkably successful in explaining the observed large 
scale structure down to length scales of ~ 1 Mpc. 

Among the most well-motivated CDM candidates 
are weakly- interacting massive particles (WIMPs), with 
masses of the order of the weak scale My/ ~ 100 GeV and 
interaction cross sections a ~ g 2 M v } 2 . WIMPs emerge 
naturally from several well-motivated particle physics 
frameworks and include the lightest supersymmetric par- 
ticle (LSP) in i?-parity conserving supersymmetry mod- 
els , the lightest Kaluza-Klein state in models with uni- 
versal extra dimensions 0, and branons in brane- world 
models 0. In addition, WIMPs naturally have thermal 
relic densities of the desired order of magnitude. 

Despite its considerable successes, however, CDM ap- 
pears to face difficulty in explaining the observed struc- 
ture on length scales < 1 Mpc. Numerical simulations 
assuming CDM predict overdense cores in galactic ha- 
los , too many dwarf galaxies in the Local Group , 
and have trouble producingenough disk galaxies without 
angular momentum loss |5j, |jj . Although there is not cur- 
rently consensus that the small scale problems of CDM 
are insurmountable the number and variety of prob- 
lems put considerable pressure on CDM and have moti- 
vated many alternative dark matter candidates. These 
include self- interacting dark matter 0|, collisional dark 
matter jjjj , thermal warm dark matter (WDM) ^lj , an- 
nihilating dark matter 0] > non-thermal WIMP produc- 
tion 0], and other proposals, such as the possibility of 
a broken scale invariance in the power spectrum |l4| . 

A common feature of these new hypotheses is that 
they preserve the successes of standard CDM on large 
scales, but reduce power on small scales. Unfortunately, 



this virtue is achieved at a cost: in contrast to WIMPs, 
these candidates are generally not well-motivated inde- 
pendently by particle physics, and their relic abundance 
is also not naturally in the correct range. For example, 
to explain the observed small scale structure, thermal 
WDM particles must have a mass greater than about 
500 eV [[IIS. 

On the other hand, the observed relic 
density is naturally achieved for masses ~ 10 eV. To 
resolve this discrepancy requires either an unreasonably 
large number (~ 10 3 ) of light degrees of freedom at the 
time of decoupling or, alternatively, a nonstandard cos- 
mology with a large injection of entropy at late times. 

Here we consider superweakly-interacting massive par- 
ticle (superWIMP) DM. In superWIMP scenarios, a 
WIMP freezes out as usual, but then decays to a sta- 
ble DM particle that interacts super-weakly [lT], ll3 • Ex- 
amples of superWIMPs include non-thermally produced 
weak-scale gravitinos [13, E E3, E3 > axinos [13], and 
quintessinos [2^ | in supersymmetry, and Kaluza-Klein 
graviton and axion states in models with universal extra 
dimensions |24|. SuperWIMPs preserve WIMP virtues: 
they exist in the same well-motivated frameworks and 
naturally have the right relic density, since they inherit 
it from late-decaying WIMPs. This latter property and 
the effect on small scale structure discussed here are ab- 
sent for thermally-produced gravitinos. 

In contrast to WIMPs, superWIMPs are produced 
with large velocities at late times. For example, gravitino 
or Kaluza-Klein graviton superWIMPs are naturally pro- 
duced at t x ~ M^JM^ - 10 5 s - 10 8 s, where the re- 
duced Planck mass M P1 = {%irG N )~ 1/2 ^ 2.4 x 10 18 GeV 
enters because these superWIMPs interact only gravita- 
tionally. This has two effects. First, the velocity dis- 
persion reduces the phase space density, smoothing out 
cusps in DM halos. Second, such particles damp the 
linear power spectrum, reducing power on small scales 
and improving consistency with structure formation. As 
we will show, these effects are sufficiently strong that 
superWIMPs may provide a natural resolution to small 
scale structure problems. Similar conclusions have been 
reached in the more general setting explored in Ref . • 
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We first consider effects coming from the veiocity dis- 
persion. Such effects may be characterized by Q = 
p/(v 2 ) 3 / 2 , the dark matter mass per unit volume of 6- 
dimensional phase space, where p is the mass density 
and {v 2 ) is the velocity dispersion [2(|. Q has a number 
of important properties. The fine-grained value of Q re- 
mains constant for collisionless, dissipationless gases. It 
may be determined analytically in scenarios with decay- 
ing dark matter. In addition, the coarse-grained value 
of Q can only decrease, a property that follows from the 
relation of Q to thermodynamic entropy. 

The coarse-grained value of Q can be estimated from 
rotation curves, gas emission, and gravitational lcns- 
mg Galaxies with coarse-grained Q near Qq = 

1.0 x 10" 27 GeV 4 ~ 1.2 x 1(T 4 (A/ Q /pc 3 )/(km/s) 3 have 
been observed |27|. Given the properties of Q noted 
above, this imposes a lower limit on fine-grained Q of 
Q > Qq. At the same time, it has been argued |2r3. l27j 
that values of Q close to Qo with, for example, Qo < Q < 
4Qoi are preferred, as they reduce the maximal cuspincss 
of galactic halos to that actually observed. 

We now determine the fine-grained value of Q in super- 
WIMP scenarios. Throughout this work, we assume that 
superWIMPs are produced between 10 3 s and 10 12 s so 
that the Universe is radiation-dominated and the number 
of effective relativistic degrees of freedom is constant. We 
define a relativistic version of Q, Q = p/{u 2 ) 3 ^ 2 , where 
u = p/m is the three- momentum normalized by the parti- 
cle's mass. Because u redshifts, Q has the advantage that 
it is invariant given expansion in the Universe for both 
relativistic and non-relativistic matter, but reduces to Q 
at late times when the matter becomes non-relativistic. 

We may estimate the value of Q at structure forma- 
tion by determining the value of Q when superWIMPs 
are produced. For simplicity, we assume that all su- 
perWIMPs are produced at the decay lifetime Tx- The 
exponential distribution of production times gives 0(1) 
corrections to the results described here and will be dis- 
cussed in detail elsewhere 25, 281. We find 
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where ux = u(tx)- Note that for ux ~ 1 and tx ~ 10 6 s, 
natural values in the cases of gravitino and Kaluza-Klein 
graviton superWIMPs, the phase space density is in the 
preferred range to eliminate cuspy halos. 

We now turn to the effect on the power spectrum. 
Initially, the matter density has small inhomogeneities. 
These inhomogeneities evolve linearly at first, but even- 
tually evolve non-linearly to form the structure observed 
today. In the present Universe, the scale of non-linearity 
is expected to be around 30 Mpc. After this point, TV- 
body simulations are necessary to analyze the evolution 
of the power spectrum. These analyses show that CDM 
predicts an excess of power on scales under 1 Mpc. 



One way to reduce power on small scales and ame- 
liorate this problem is through the free-streaming of 
DM. From its production time tx until matter-radiation 
equality at rj E Q ^ 2.2 x 10 12 s, superWIMPs can stream 
out of overdense regions into underdense regions, smooth- 
ing out inhomogeneities. A free-streaming scale much 
larger than 1 Mpc is excluded by observations of Ly- 
man alpha clouds jl^- (For instance, for WDM, the con- 
straints mwDM ^ 550 eV ^| and towdm }t 750 eV 
correspond, given our definition of free-streaming scale, 
to Aps S 1-4 Mpc and Aps < 1.0 Mpc, respectively.) 
However, values close to this (roughly, 1.0 Mpc > Afs ^ 
0.4 Mpc) could resolve the present disagreements. 

The free-streaming scale for superWIMP dark matter 
can be estimated to be 
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where a(t) is the cosmic scale factor, v(t) is the super- 
WIMP velocity, and 
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where ffcp = ct E ct/a(tEQ) — 93 Mpc and ueq = 
w(*eq) 03113 In the common case that superWIMPs 
are relativistic when produced (u x ^ 1) but non- 
relativistic at iEQ (meq < 1)) Eq. simplifies to 
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demonstrating that production times of ~ 10 6 s naturally 
also provide the preferred order of magnitude for the free- 
streaming scale. Note that in this case, Q and Afs are 

1/2 

both functions of uxTx , and so are correlated, as they 
are for WDM. However, in the superWIMP scenario with 
ux ^ 1) Afs < A(t^Q) and this degeneracy is broken, 
opening up new possibilities [28j |. 

We now consider the particular case of gravitino su- 
perWIMPs. The next-to-lightest supersymmetric parti- 
cle (NLSP) may be either a scalar (e.g., a sneutrino or 
a slepton) or a fermion (a neutralino). If the NLSP is a 
sneutrino, gravitinos are produced at time |18| 
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with ux = {mf, — m 2 =,)/(2777 J >?77 ( =,). These then deter- 
mine Q and Afs through Eqs. CI - ©; the results are 
shown in Fig. ^ Remarkably, preferred values of Q and 
Afs are simultaneously realized in superWIMP scenarios 
with natural, weak-scale masses with ttinlsp ^ 500 GeV. 
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FIG. 1: Preferred regions (shaded) of phase space density Q 
and free-streaming length Afs in the (mswiMP, Am) plane, 
where Am = m^LSP — ^swimp, for gravitino superWIMPs 
with a sneutrino NLSP. The regions under both bands are 
disfavored. In the regions above both bands, superWIMP 
dark matter becomes similar to CDM; representative values 
of Q and Afs are shown. Contours of typical lifetimes re are 
also shown. We have assumed fiswiMpft 2 = 0.11. 



The superWIMP scenario is also constrained by Big Bang 
nucleosynthesis and the Planckian spectrum of the cos- 
mic microwave background. These constrain ts [301 have 
been evaluated in several studies 0, 0, 0, l2dL l2lL l3l| . 
For a v NLSP, the main constraint comes from hadronic 
and electromagnetic energy produced in the three-body 
decays v — > IWG, vZG. These modes have small branch- 
ing ratios of the order of 10~ 3 and so the parameter space 
is not strongly constrained for tx ^ 10 6 s [2(3 ]. 

For a charged slepton NLSP, the lifetime is identical 
to that given in Eq. JSJ- However, charged sleptons 
are not collisionlcss and are electromagnetically coupled 
to the baryon-photon plasma. The opposite tendencies 
of pressure repulsion and gravitational attraction gen- 
erate acoustic waves with density perturbation oscilla- 
tions of photons, baryons and sleptons. After gravitino 
production, however, the photon-baryon fluid is coupled 
only gravitationally to the neutral superWIMP. Power is 
therefore reduced on scales that enter the horizon before 
this decoupling, and this effect can be more important 
than the free-streaming damping discussed above |28ll32| . 

If the decaying WIMP is a neutralino, the superWIMP 
production time is different. In general, the neutralino 
is a mixture of neutral Bino, Wino, and Higgsino states. 
For the specific case of a photino, the lifetime is 0] 
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FIG. 2: Same as in Fig. Q but for gravitino superWIMPs 
with a photino NLSP. 



The resulting Q- and A^s-preferred regions are shown in 
Fig. As in the sneutrino NLSP case, preferred val- 
ues for both quantities are obtained for natural weak- 
scale masses; if anything, the preferred superWIMP and 
NLSP masses are lighter and more natural. For typical 
neutralinos, hadrons are produced in two-body decays, 
leading to extremely severe constraints from Big Bang 
nucleosynthesis. For photinos, however, the constraints 
are relaxed, since hadrons are produced only in three- 
body decays, and the superWIMP resolution to small 
scale structure problems may be realized. 

In conclusion, we have examined the implications of su- 
perWIMP dark matter for small scale structure. Because 
superWIMP dark matter is produced with large veloc- 
ity in late decays, its phase space density is decreased, 
smoothing out halo cusps. At the same time, super- 
WIMPs damp the linear power spectrum, which lowers 
the concentration of dark matter in the cores of galactic 
halos. This effect may bring numerical simulations into 
agreement with observations of dwarf galaxies and reduce 
the excess of dwarf galaxies relative to CDM predictions. 

The effects on small scale structure may also be 
achieved by WDM or dark matter with exotic interac- 
tions. In contrast to those possibilities, however, su- 
perWIMPs are automatically present in particle physics 
models with supersymmetry or extra dimensions and are 
naturally produced with the correct relic density. Dark 
matter produced in late decays will necessarily be warmer 
than CDM. It is remarkable, however, that for super- 
WIMP gravitinos, where the production times and ve- 
locities are determined by the fixed energy scales Mew 
and Mpi, the predicted values of both Q and Afs ar e in 
favorable ranges to resolve outstanding problems with- 
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out violating other constraints from cosmology and par- 
ticle physics. SuperWIMPs therefore appears to com- 
bine the most appealing features of both CDM and 
WDM. This explanation will be probed by future ob- 
servations, especially those constraining the epoch of 
reionization |25l l33j| : reionization by redshift 6 implies 
Q > 0.1 x Qo, compatible with the preferred values ana- 
lyzed in this work, but confirmation of indications from 
WMAP of earlier reionization could restrict parameter 
space greatly. At the same time, given the virtues de- 
scribed here, it would be especially interesting to see if 
the promise of super WIMPs is realized by Y-body simu- 
lations of structure formation and semi-analytic analyses. 
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